Dissipation and elliptic flow at RHIC 
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We compare elliptic flow evolution from ideal hydrodynamics and covariant parton transport 
theory, and show that, for conditions expected at RHIC, dissipation significantly reduces elliptic 
flow even for extreme parton cross sections and/or densities a gg x dN / dr/(b—0) ~ 45 mb x 1000. 
The difference between transport and hydrodynamic elliptic flow is established rather early during 
the evolution of the system, but the buildup of elliptic flow is surprisingly insensitive to the choice 
of the initial (formation or thermalization) time in both models. 
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Introduction. During the early stage of ultra- 
relativistic heavy ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) and the future Large Hadron Col- 
lider (LHC), a deconfined phase called the quark-gluon 
plasma (QGP) is expected to be formed. One way to in- 
fer the properties of this dense partonic phase is to study 
its collective behavior. 

An important experimental probe of collective dy- 
namics in noncentral A + A reactions is elliptic flow, 
t>2 = (cos(2(/>)), the second Fourier moment of the az- 
imuthal momentum distribution Recent data from 
RHIC for Au+Au at i/IWn ~ 130 - 200 GeV show a 
large anisotropy of particle production in the transverse 
plane with V2 values reaching up to 0.2 and saturatingin 
the transverse momentum region 2 < p± < 6 GeV |2J, y( . 

For p± < 2 GeV, i.e., for the bulk of the particles pro- 
duced at RHIC, the elliptic flow data can be reproduced 
from ideal (Euler) hydrodynamics. In this completely 
nondissipative theory elliptic flow is sensitive to the equa- 
tion of state (EOS) of dense nuclear matter 0,0>IEI3- 
It is remarkable that agreement with the data requires 
HI IE an EOS with deconfinement phase transition, 
providing one of the strongest arguments for QGP for- 
mation at RHIC. 

Despite its successes, the assumption of no dissipa- 
tion has limitations. For example, ideal hydrodynam- 
ics fails || to saturate elliptic flow for p± > 2 GeV, 
and overpredicts [2( the "long" and "out" pion Hanbury- 
Brown and Twiss (HBT) interferometry radii. Both of 
these shortcomings are likely due to the neglect of dis- 
sipation 0, 0, |3 E3 caused by nonzero mean free 
paths. In any case, the mere proof that the effect of 
dissipation is negligible (for some observable) already re- 
quires a framework that allows for dissipation. 

For systems near equilibrium, dissipative corrections 
can be studied via Navier-Stokes (viscous) hydrodynam- 
ics. In this theory, the evolution is determined by the 



equation of state, the shear and bulk viscosities, and 
the heat conductivity. Even though a general Lorentz- 
covariant formulation exists, only 1+1D relativistic so- 
lutions are known 0, [^. Nevertheless, recent esti- 
mates [ill E2I Il5j | find considerable viscous corrections 
to the evolution of the system, and to spectra, elliptic 
flow V2(p±), and pion HBT parameters at RHIC. 

To study dissipative effects arbitrarily far from equilib- 
rium, one can utilize covariant parton transport theory 
El H El El El El In this approach elliptic flow 
depends mainly on the effective scattering cross section of 
partons produced in the collision jlOl [lfj . Both the mag- 
nitude and the saturation feature of elliptic flow at RHIC 
can be reproduced E3 if very large (~ 45 mb) elastic 
parton-parton cross sections are assumed. It is puzzling 
that these values are an order of magnitude above con- 
ventional perturbative QCD (pQCD) estimates. There 
are, however, recent proposals to alleviate this prob- 
lem mim. 

A popular interpretation of the joint success of ideal 
hydrodynamics and parton kinetic theory at RHIC is that 
with ~ 45 mb elastic parton-parton cross sections the 
opacities are large enough to reach the ideal hydrody- 
namic limit, at least at px < 2 GeV. However, that con- 
clusion ignores that calculations based on the two theo- 
ries corresponded to different initial conditions and ther- 
modynamic properties. 

In this paper we provide a systematic comparison of 
ideal hydrodynamics and transport theory and investi- 
gate the importance of dissipative effects at RHIC. Be- 
cause a thermodynamically consistent microscopic de- 
scription of the QGP phase transition is still a difficult 
open problem, we gain insight via studying an ideal gas 
of massless gluons. 

Precursors of this study in Refs. [pi fisj l analyzed the 
transverse energy of particles and found ~ 20% dissipa- 
tive corrections. Here we focus on differential elliptic flow 
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and show that it is more sensitive to dissipation. 

Ideal hydrodynamics and covariant transport theory. 
Ideal (Euler) hydrodynamics is a convenient Lorentz- 
covariant dynamical theory formulated in terms of 
macroscopic quantities. In the context of heavy ion col- 
lisions, its main advantage is the ability to treat phase 
transitions, while its main limitation is that it is only 
applicable to systems in local kinetic equilibrium. 

The hydrodynamical equations of motion are the local 
conservation laws of energy-momentum and net charge 

d^{x) = 0, d M N£(x) = (1) 

[x = (t,x) is the Minkowski four-coordinate]. In ideal 
hydrodynamics the energy-momentum tensor is assumed 
to be that of an ideal fluid, T^ v = (e + p)u^u v - pg*", 
where e(x), p(x), and are the local energy density, 

pressure, and flow velocity. The charge current is related 
to the local charge density n c via N£ — n c u^ . Once 
the equation of state p(e, n c ) is specified, Eqs. (JIJ can be 
utilized to follow the evolution of the system from any 
given initial state in local kinetic equilibrium. 

For heavy-ion collision applications, hydrodynamics 
has to be supplemented with a freezeout description be- 
cause the assumption of local equilibrium breaks down 
as the local mean free path A = \/[a{s)n(xj\ becomes 
comparable to the Hubble radius during the expansion. 
A common approach is to assume that microscopic scat- 
tering rates T sc (T,n) ~ o~(T)n drop so quickly that one 
can consider a sudden transition from local equilibrium 
to a noninteracting gas. In this case, the fluid is con- 
verted to particles on a 3D spacetime hypersurface via 
the Cooper-Frye formula [22^. For systems in chemical 
equilibrium, a constant temperature (or constant energy 
density) hypersurface is usually chosen because the only 
relevant parameter is the temperature, T SC (T, n(T)). 

Covariant parton transport theory, in contrast to hy- 
drodynamics, is based on microscopic quantities, namely, 
the phasespace distributions of (quasi)particles and mi- 
croscopic transition probabilities. The main advantage 
of this approach is that it is applicable out of equilib- 
rium and models freezeout self-consistently. However, it 
cannot describe phase transitions (without coupling to 
classical fields). 

We consider here, as in Refs. [3 III III III III E3, the 
simplest but nonlinear form of Lorentz-covariant Boltz- 
mann transport theory in which the on-shell phase space 
density f(x,p), evolves with an elastic 2^2 rate as 

tfdpfi = S(x,ft) + Y^JJJ (/a/4 - hh) \Mi2^i\ 2 

234 

x 5 i {p l +p 2 -P3-Pi) -(2) 

Here |.M| 2 is the polarization averaged scattering ma- 
trix element squared, the integrals are shorthands for 
^ = J d 3 Pi/(2Ei), while f 3 = f(x,p ). The initial con- 
ditions are specified by the source function S. For our 



applications below, we interpret f(x,p) as describing an 
ultrarelativistic massless gluon gas with g = 16 degrees 
of freedom (8 colors, 2 helicities). 

Eq. (J2J, which corresponds to Boltzmann statistics, 
could in principle be extended for bosons and/or for in- 
elastic processes, such as gg <-> ggg. However, no prac- 
tical covariant algorithm yet exists that can handle, at 
the opacities expected at RHIC, the new nonlinearities 
these extensions introduce. We therefore limit our study 
to quadratic dependence of the collision integral on /. 

Boltzmann's H -theorem states that Eq. © drives the 
system towards a fixed point, global equilibrium. In the 
hydrodynamic limit, i.e., when \A4\ 2 — * oo, the transport 
evolution approaches the ideal hydrodynamic evolution, 
f(x,p) = gexp[(fi{x) - p v u v (x))/T{x)]/{2Tif, provided 
hydrodynamics is not frozen out via some prescription. 
However, for a finite matrix element, an expanding sys- 
tem, if it ever equilibrates, sooner or later evolves out of 
equilibrium (except for very special examples |2fljV 

Numerical results. We study elliptic flow in a typical 
midperipheral Au+Au collision at RHIC with impact pa- 
rameter 6 = 8 fm. The initial conditions were taken from 
|lfj|: a longitudinally boost invariant Bjorken tube in lo- 
cal kinetic equilibrium at temperature To at proper time 
To fm/c, with a transverse gluon density distribution that 
is proportional to the binary collision distribution for two 
Woods-Saxon nuclei. A thermal momentum distribution 
is necessary because hydrodynamics is limited to such 
initial conditions. Unlike in usual hydrodynamic calcu- 
lations, the initial temperature was the same in the en- 
tire system even though the density of particles was not. 
This means that the system was not in chemical equilib- 
rium. Motivated by parton-hadron duality |24| , the initial 
gluon density was normalized to dN g /di](b = 0) = 1000 
to match the observed dN c h/dr] ~ 600. Accordingly, 
dN g /dri w 250 for b = 8 fm. 

We solved Eq. (01 via the MPC algorithm^, which 
utilizes parton subdivision^! Il8l | to maintain Lorentz 
covariance. The elastic parton cross section was isotropic 
with o~t ot = 3, 8, and 20 mb, which via transport opac- 
ity scaling |10j is approximately equivalent to using a 
Debye-screened leading-order pQCD gg — > gg cross sec- 
tion da/dt cx l/(t — ^ 2 ) 2 with a tot ~ 7, 19, and 47 mb. 
Parton subdivisions of 180 were used to stabilize the nu- 
merical results. 

The hydrodynamic solutions were obtained via the 
code used in Refs. 6, 26] . Particle number was conserved 
explicitly because we are comparing to transport theory 
with elastic 2^2 interactions. For the massless par- 
ton gas considered here, the EOS is e = 3p. Because the 
cross section is energy-independent, freezeout was chosen 
to occur at a proper density of n g — 0.365 fm~ 3 . This 
corresponds to the gluon density of a chemically equili- 
brated gluon gas at T = 120 MeV. 

Figure H shows elliptic flow V2{p±) at freezeout for 
Tq = 700 MeV and tq = 0.1 fm/c (which are the param- 
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FIG. 1: Comparison of V2(p±) from ideal hydrodynamics 
(solid curve) and transport theory as a function of the parton 
cross section (dashed curves), for tq — 0.1 fm/c. 
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FIG. 2: Time evolution of V2(p±) from ideal hydrodynamics 
(solid curves) and transport theory for a gg = 47 mb (dashed 
curves) . 



eters used in jlOj ) . From identical initial conditions and 
thermodynamic properties (e = 3p EOS) ideal hydrody- 
namics produces a much larger anisotropy than transport 
theory for all three parton cross sections a gg = 7, 19, and 
47 mb. Even for the largest ~ 45 mb parton cross section 
studied, dissipation reduces vi by 30 — 40%. This effect is 
larger by a factor 1.5-2 than what earlier studies [fIUs! 
found for the final transverse energy dEr/dy |27j . 

The reason that both ideal hydrodynamics and par- 
ton transport could reproduce the same RHIC vi data 
is that the earlier hydrodynamical calculations utilized a 
softer EOS (e > 3p) that incorporates the QGP phase 
transition, which reduces elliptic flow. Though dissipa- 
tive effects for such realistic equations of state are not 
yet calculable, there is no indication that they would be 
smaller than for the ideal gas studied here. 

The difference between transport and hydrodynamic 
elliptic flow is established rather early during the evolu- 
tion, as shown in Fig. where we plot V2(p±) computed 
over hypersurfaces of constant r. Already by r w 1 fm/c, 
the hydrodynamic anisotropy is a factor two or more 
above the transport, especially at high p±. The transport 
V2 builds up smoothly and at much the same pace both 
at high and low p±, and by r = 3 fm/c most (~ 80%) of 
the final anisotropy is there. This reinforces the short 
few- fm/c timescales found for the pj_-integrated v 2 in 
0, 0| . On the other hand the hydrodynamical develop- 
ment of anisotropy is more pj_-dependent. Above pj_ = 1 
GeV most (85-90%) of the anisotropy is also developed 
at r = 3 fm/c, but at low p±, e.g. at pj_ = 0.5 GeV, 
anisotropy still increases by a factor of 1.5 before freeze- 
out at r = 5.3 fm/c. 

A qualitative explanation for the remarkable initial lag 
in the vq. evolution for the transport is that for an equi- 
librium initial condition Eq. PJ in fact corresponds to 
free streaming because the collision term vanishes ex- 



actly. Momentum observables can start changing only 
after the system streams out of equilibrium. In contrast, 
in hydrodynamics, the initial pressure gradients induce 
changes in the flow pattern and therefore immediately 
start to generate elliptic flow. 

Because of scalings of the dynamical equations, the 
above results are quite general. For the transport, the so- 
lutions depend only on three scales^3 : To, a trivial scale 
that fixes the momentum units; and two nontrivial scales 
atotdN/dr] and R/tq. For ideal hydrodynamics, one can 
similarly prove that three scales apply: To (momentum 
units), dN/drj (linear scale for particle number, e.g., the 
spectra), and R/tq. For example, for twice larger To and 
one-third the dN/dr], the transport V2(p±) solutions are 
the same provided Otot is increased three times and the 
p± axis is stretched by factor of two, while the hydrody- 
namic V2(pi_) solutions can be obtained via the stretching 
of p± axis alone if the freeze-out density is scaled in the 
same way as the initial density. 

Another important difference between the earlier hy- 
drodynamic and transport calculations was that the hy- 
drodynamic evolution started at an assumed thermaliza- 
tion time tq — 0.6 fm/c, whereas the transport evolution 
started at a formation time estimate of tq = 0.1 fm/c. 
One might naively expect that the much smaller initial 
time for the transport leads to more rapid departure from 
equilibrium and stronger dissipative effects. However, 
for both theories elliptic flow v 2 (p±) is largely insensitive 
to the initialization time, as shown in Fig. [3] where we 
compare the above mentioned results for to = 0.1 fm/c 
and T = 700 MeV with those for t = 0.6 fm/c and 
To = 385 MeV. To approximately preserve typical par- 
ticle momenta at freezeout, we have rescaled the initial 
temperature according to a ID Bjorken expansion, under 
which T - r" 1 / 3 . 

The main reason for the insensitivity to initial time 



4 




P± [GeV] 

FIG. 3: Weak initialization time dependence of V2(p±.) from 
ideal hydrodynamics (pluses) and transport theory for a gg — 
47 mb (boxes). 

is that the only relevant parameter that differs between 
the two calculations, tq/R, is in both cases much smaller 
than one. At early times when (r — tq) <C R, the sys- 
tem essentially expands only longitudinally. Therefore, 
between r = 0.1 and 0.6 fm/c very little elliptic flow is 
generated, as can be seen in Fig. |21 Our results reinforce 
similar findings for ideal hydrodynamics by Ref. Q, and 
also show that the elliptic flow pattern is largely insensi- 
tive to the initial time even in the presence of significant 
dissipation, as long as tq <C R. 

Conclusions. In this paper we studied elliptic flow 
V2(j>±) using ideal hydrodynamics and parton transport 
theory for an ideal gas of massless partons. From identi- 
cal initial conditions and thermodynamic properties, dis- 
sipation significantly reduced V2(p±), even for extreme 
~ 45 mb elastic parton cross sections that are of the or- 
der of hadronic cross sections. In addition, comparison 
to earlier works shows that elliptic flow is more sensitive 
to dissipation than the final transverse energy. 

These results indicate that the large elliptic flow seen 
at RHIC can be interpreted in (at least) two ways: either 
as i) a harder nuclear equation of state and strongly dissi- 
pative evolution, or ii) a softer equation of state and neg- 
ligible dissipation which means that the system stays in 
essentially perfect local kinetic equilibrium via a mecha- 
nism whose microscopic origins we do not yet understand. 

We emphasize that our results correspond to the simple 
ideal gas equation of state. A detailed investigation of 
dissipation in nuclear collisions will have to consider in 
the future a more realistic nuclear equation of state that 
includes the hadronization phase transition. 

Finally, though our transport solutions do not reach 
the ideal hydrodynamic limit, it would be very interesting 
to check how far they are from the Navier-Stokes limit. 
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